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ABSTRACT 

A two component model of nonthermal dark matter is formulated to simultaneously explain the 
Fermi-LAT results indicating a 7-ray excess observed from our Galactic Centre in the 1-3 GeV 
energy range and the detection of an X-ray line at 3.55 keV from extragalactic sources. Two 
additional Standard Model singlet scalar helds S'2 and S'3 are introduced. These helds couple 
among themselves and with the Standard Model Higgs doublet H. The interaction terms among 
the scalar helds, namely H, S 2 and S 3 , are constrained by the application of a discrete Z2 x Z2 
symmetry which breaks softly to a remnant II 2 symmetry. This residual discrete symmetry is 
then spontaneously broken through an MeV order vacuum expectation value u of the singlet 
scalar held S 3 . The resultant physical scalar spectrum has the Standard Model like Higgs as 
with ~ 125 GeV, a moderately heavy scalar with 50 GeV < My.^ < 80 GeV and a 
light with ~ 7 keV. There is only tiny mixing between y^ and y^ as well as between 
y^ and y^. The lack of importance of domain wall formation in the present scenario from the 
spontaneous breaking of the discrete symmetry provided u < 10 MeV, is pointed out. We 
hnd that our proposed two component dark matter model is able to explain successfully both 
the above mentioned phenomena — the Fermi-LAT observed 7-ray excess (from the y^ bb 
decay mode) and the observation of the X-ray line (from the decay channel y^ —)■ 77) by the 
XMM-Newton observatory. 
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1 Introduction 


The presence of Dark Matter (DM) in the Universe is now an accepted reality. So far, its 
existence has been inferred only from its gravitational effects. The latter include rotation curves 
of galaxies, gravitational lensing, observations of the Bullet cluster etc. However, a very attractive 
proposition is that DM consists of Weakly Interacting Massive Particles (WIMPs) [1] and this 
is the paradigm that we adopt. The possibility of detecting DM WIMPs through direct and 
indirect processes is now being vigorously pursued by experimenters. Meanwhile, the amount 
of DM present in the Universe has been precisely determined by means of results from the 
PLANCK satellite [2] whose instruments have probed and analysed anisotropies in the smooth 
cosmic microwave background. One can think of two possible scenarios for the production of DM 
particles in the early Universe. (1) This could have been due to thermal processes with Standard 
Model particles interacting with the thermal plasma in the expanding soup ball whereby all 
possible particle-antiparticle pairs were produced in a reverse chemical process. The decoupling 
of the DM particles occurred when their interaction rates fell short of the expansion rate of the 
Universe. Being out of equilibrium, the dark matter particles “froze” to a particular relic density. 
(2) The DM particles could have originated nonthermally. In this mode, they got produced from 
out of equilibrium decays of heavier particles. That could have occurred at the stage of the 
preheating of the Universe from the inflaton energy. 

Dark Matter, if accumulated in considerable measure at highly dense regions of celestial bod¬ 
ies, such as at the Galactic Centre (GC) or other sites in a galaxy, may undergo self annihilation. 
That would lead to high energy photons or fermion-antifermion pairs. Such hard photons may 
also arise from the radiative decay of a metastable DM particle such as a sterile neutrino or a 
scalar. Therefore, the observation and analysis of such high energy photons could not only lead 
to an indirect detection of DM particles but also provide insights into the physics of Dark Matter. 

Recently, there have been two major observations of such high energy photons: y-rays from 
the centre of our Milky Way galaxy and X-rays from other galaxies and galaxy clusters. The 
observation of a weak unidentified line in the X-ray spectrum obtained from the XMM-Newton 
observatory, has been reported by Bulbul et al [3], from an analysis of data from 73 galaxy 
clusters including Perseus and others. It is also reported in same article [3] that the best fit 
value of the observed X-ray flux obtained from the XMM-Newton MOS(PN) observation is 
X 10“®(3.9l5;o ^ 10“®) photons cm“^ s“^ while the energy of the X-ray line is E = 
3.57 ± 0.02(3.51 ± 0.03) keV. This has been conhrmed later by Boyarsky et al [1]. In Ref. |1] 
this 3.55 keV line has been claimed to have been observed from the Andromeda as well as other 
galaxies in the Local Group. Separately, analyses of the Fermi-LAT data by several groups 0-II1 
during the last few years have been indicating the presence of a significant excess at an energy 
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range 1-3 GeV in the 7 -ray spectrum observed from regions close to the centre of our Milky 
way galaxy. There have been attempts to explain the 3.55 keV line from other astrophysical 
processes. For example, in Ref. [12] the authors claimed from their analysis of XMM-Newton 
data that a 3.55 keV line from the Galactic Gentre could be interpreted from a known plasma 
line of astrophysical origin. In another work by Garlson et. ah [16] the 3.55 keV line from the 
Perseus cluster is related to the cool core of the cluster. There have also been attempts m 
to obtain the XMM-Newton X-ray emission line from dwarf galaxies. Authors have considered 
sterile neutrino dark matter to explain the signature of this emission and accordingly have given 
upper bounds on the mixing angle of the relevant sterile neutrino DM. X-ray emission lines of the 
order of 3.5 keV have also been probed from galaxy clusters, the Galactic Gentre and M31 where 
analyses [18] have been made considering sterile neutrino dark matter with a few keV mass. But 
the authors could not report any ~ 3.5 keV excess. For the case of the 1-3 GeV 7 -ray excess 
from the Galactic Gentre, the authors of Ref. [TU] attempted to demonstrate the origin of the 
spectral and the angular features of this excess from inverse compton scattering of high energy 
electrons available from a burst event in the distant past. A possible millisecond pulsar origin 
of this 7 -ray excess from the Galactic Gentre has been discussed in Ref. [2U]. But none of the 
above interpretations has addressed both the phenomena of the X-ray line and the 7 -ray excess 
together. Hence it is natural to conjecture a common origin of both phenomena from processes 
involving DM particles. There have been many dark matter models in the literature explaining 
either the Fermi-LAT observed 7 -ray excess [21]-[13] or the X-ray line [12]-[68] detected by the 
XMM-Newton observatory. However, we are among the hrst two [22] to propose a new single 
dark matter model explaining both the phenomena simultaneously, although our model of two 
SM-singlet scalars is very different from that of Ref. [69]. While our explanation of the 7 -ray 
excess is new, our explanation of the X-ray line uses the same mechanism as that of Ref. 

Once again, there have been previous works on two component dark matter models mi, |7D1-|71, 
however none of those has been aimed at explaining the Fermi-LAT 7 -ray excess and the 3.55 
keV X-ray line simultaneously like we do. 

In this work we propose a single two-component model for nonthermal dark matter applying 
it to two completely different parts of the electromagnetic spectrum. We have two distinct scalar 
dark matter particles, one being light and hence “warm” (~ 7 keV in mass) and another which 
is moderately heavy and “cold”, being in the mass range 50—80 GeV. Thus, we extend the 
scalar sector of the Standard Model (SM) of particle physics by two real scalar field S 2 ans As, 
both of which are singlets under the Standard Model gauge group SU(2 )l x U(1)y- We envisage 
nonthermal WIMPs and propose that both the DM particles were produced from decays of the 
Standard Model Higgs boson as well as from the pair annihilation of SM particles such as fermions, 
gauge bosons and Higgs bosons. The production processes of these DM particles took place in 
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the early stage of the Universe when its temperature fell below the electroweak phase transition 
scale. We need not consider the chemical equilibrium in the thermal bath before electroweak 
symmetry breaking but are concerned with the nonthermal processes occurring afterwards that 
lead to the creation of the two lighter scalars. The interactions of the new scalar helds S 2 , S 3 — 
among themselves as well as with the SM Higgs field H — are controlled by appropriate discrete 
symmetries which are softly broken down to a residual Z 2 symmetry under which both S 2 and 
S 3 are odd while all other fields are even. This residual discrete symmetry gets spontaneously 
broken when S 3 develops a Vacuum Expectation Value (VEV) u of order MeV. This results in a 
3x3 mass squared matrix in the three dimensional space of the residual neutral SM Higgs held 
and the two new singlet scalar helds. The diagonalisation of that mass squared matrix leads to 
three physical scalar masses which are taken as ~ 125 GeV, M-y ~ 50-80 GeV and My ~ 

f J Xi ’ X2 X3 |-| 

7 keV corresponding respectively to the SM-like Higgs the two DM particles and x^ 

The comoving number densities and of at the present epoch are calculated 

by numerically solving the corresponding two coupled Boltzmann equations describing their 
temperature evolution. These equations take into the account the roles played by the decays 
y^ —)■ XjX^ (j = 2, 3) as well as the pair annihilation processes xx X-X- where x can be 
Z, /(/), y^, y^,/being any SM fermion. The individual relic densities and 

at the present temperature then follow in a straightforward way once the above are computed. 
The temperature variation of these, as well as of the total relic density are numerically 

studied for three cases: > Dv , and ~ . The dependence on the mass 

, which is not pinpointed, is also studied. The total relic density is always found to lie in 
the range 0.1172 < St'YhS < 0.1226 in conformity with the latest PLANGK results [2]. 

We further consider the issue of domain wall formation from the restoration of the discrete 
symmetry Z 2 and argue, a la Babu and Mahapatra IS], to conclude that the corresponding 
energy density would be too little to have any effect on the possible overdosing of the Universe 
or the near-isotropy of the GMB so long as u is bounded from above by ~10 MeV. Finally, we 
take up the two main issues towards which this paper is addressed: the 7 -ray excess observed 
by Fermi-LAT at 1-3 GeV energies from our GG and the detection of an anomalous 3.55 keV 
X-ray line from extragalactic sources by the XMM Newton observatory. For the former, we argue 
that y^ —)■ bb is the dominant decay mode acting as a source of GeV energy 7 -rays from decays 
of neutral pions created during the hadronisation of the b and the b. The computation of the 
resultant 7 -ray flux from the decay of the dark matter candidate y^ is carried through using 

'^This mass range for the dark matter component is required to explain the Fermi-LAT excess 7 -rays from 
the decay channel y^ —>■ bb (c/. Section [5] for a more detailed discussion) while the mass of the other DM 
component y^ needs to be 7.1 keV so that its decay y^ —77 produces two monoenergetic photons each with 
energy 3.55 keV. 
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the NFW [75] halo prohle near the GC. A comparison with the Fermi-LAT data provides an 

excellent £t for the mass range of Xj chosen by us. The observed 3.55 keV X-ray line arises in 

our model from the decay Xg —t 77 . We hrst compute the partial decay width of the DM 

particle x utilising three operative one loop diagrams. This partial decay width of x for the 

Ox 

channel Xg ^ 77 , however, needs to be modihed by the factor for the purpose of computing 
the differential X-ray flux since we are working in a multicomponent dark matter scenario. We 
hnd that the modified decay width for the channel Xg 77 lies within the range predicted in 
Refs. [21 m ESI so long as the VEV u of the singlet scalar held S 3 is bounded from below by 2.4 
MeV. The decay modes Xg bb and Xg 77 arise respectively through the nonzero mixing 
of both the dark matter candidates Xg ^rid Xg with the SM-like Higgs boson x^ • Moreover, the 
couplings are sufficiently small (owing to the nonthermal origin of X 2 ^md 
Xg) to evade all the existing constraints from dark matter detection experiments [7S1E7]. 

The rest of the paper is organised as follows. In Section [2] we describe the model and discuss 
its interscalar interactions as well as the emergent 3x3 mass squared matrix in the scalar sector. 
Section [3] presents the calculation of the relic densities at the present temperature 

and a discussion of their variation with temperature. In Section 0] we discuss the issue of possible 
domain wall formation and point out why it is unimportant in the present context. The spin 
independent elastic scattering cross sections of both the dark matter particles are computed in 
Section (5] Section El contains our calculation of the 7 -ray hux arising from the decay of the 
dark matter component Xg and its comparison with the available Fermi-LAT data. In Section [7| 
we compute the decay width for the channel Xg 77 which is required to produce the 

observed 3.55 keV X-ray line from the decay of the dark matter particle Xg- The hnal Section [ 8 ] 
summarises our conclusions. The two appendices flA.ll and lA.2p contain the algebraic expressions 
for the couplings and masses of the physical scalars x^, Xg ^3 • 


2 The Model 


We start with the Standard Model fields including, in particular, the SU(2 )l Higgs doublet 


H 


h+ 


( 1 ) 


\ V2 / 

Here is the residual neutral SM Higgs held, h~^ and are unphysical charged and neutral 
pseudoscalar Higgs helds while v is the Vacuum Expectation Value ~ 246 GeV. Thus 


{H) = 



( 2 ) 
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Our DM sector consists of two real scalar fields S 2 and S 3 which are singlets under the SM gauge 
group SU(2 )l X U(1)y- The Lagrangian of the spin zero sector of the present model can then 
be written in terms of the SU(2 )l x U(1)y gauge covariant derivative and a scalar potential 
function V{H, 52 , S 3 ): 


C = {D^H)\D^^H) + ]^d^S 2 d^S 2 + \d^S 3 d^S 3 - V{H, S^, S 3 ) . 


(3) 


The stability of the newly added scalars ^ 2,3 is ensured by the postulated invariance of the model 
under the discrete symmetry Z 2 x Z 2 with respect to which only these two helds transform 
nontrivially. The Z 2 x Z 2 charges of S 2 and S 3 are (-1, 1) and (1, -1) respectively, whereas the 
corresponding charge of every other held in the model is (1, 1). Thns, while S 2 {S 3 ) is odd (even) 
nnder Z 2 , the reverse is the case for S 2 {S 3 ) with respect to Z 2 . The fnnction V{H, S 2 , S 3 ) of 
Eq. (j3]) contains all possible Z 2 x Z 2 invariant interaction terms among the Standard Model 
Higgs donblet held H and the singlet helds S 2 , S 3 inclnding all allowed self interactions and mass 
terms. We write it explicitly in Eq. (|1]): 


V{H,S 2 ,S 3 ) = +'Ss^+'S(^S: 




+ A 


13 


(h^h 


V 

'2 
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+ X 32 {H^H)Si + A23525I 


u 


(4) 


The reqnired invariance nnder the symmetry Z 2 x exclndes terms snch as {WH)S 2 S 3 , S 2 S 3 , 
525 |. We also assnme that, dne to high scale physics, a ‘soft’ term 

E' = aS2S3 (5) 

gets added to V, explicitly breaking the Z 2 x Z 2 invariance down to Z 2 nnder which both S 2 and 
S 3 are odd. We have introdnced six ‘hard’ conplings in E — 2 , 3 ) -^ 12 ) A 23 , A 13 — all of which 

describe qnartic held interactions in their leading terms. Besides the six quartic conplings, we 
have also a dimensional conpling p^. The nonleading terms in the potential V have been chosen 
in a way that V is manifestly minimised at {H) = {S 2 ) = 0, {S 3 ) = u. The VEV u needs to 

be in the range 2 MeV < u < 10 MeV (see Sections 01 [7] for details). As a resnlt, the residnal 


y broken. The stability of the potential V in 
of Ref. [TS]. The reqnired conditions are: 

Eq. (01) can 

^1.2,3 ^ ^ ; 

( 6 ) 

P 2 >0 , 

(7) 

Ai2 + \/^ i ^2 ^ 0 5 

( 8 ) 

Ai3 + \/^ i ^3 ^ ^ ’ 

(9) 


( 10 ) 
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( 11 ) 


-^2 (Ai 2 + \/«i«2) (-^13 + ( 2 A 23 + \/^2^3) 

+ \/^ 1 ^ 2 ^3 + + Ais-y/^ + 2A23\/^ > 0 . 


In analogy with the residual neutral SM Higgs held s^, we can dehne the corresponding 
Sj = S 2 and = S 3 — u. Now, in the s^-s^-s^ system, the squared mass matrix is given by 




^ 0 2Ai3Mn ^ 

0 P 2 + Ai 2 n^ + ‘2\23U^ a 
y 2Ai3Mn a 2 k^v? y 


( 12 ) 


The eigenvalues of are designated and . The eigenstates, namely the physical 

ACf /A,2 ^^3 

scalars are linearly related to Si, S 2 , S 3 by an orthogonal transformation which may 

be characterised in the usual way by the mixing angles 612 , 623 and ^ 13 , The latter are the rotation 
angles needed to rotate the S'i-S' 2 , S 2 -S 3 and S' 3 -S'i mass squared submatrices ( 2 x 2 matrices) 
sequentially such that the whole mass squared matrix (Eq. flT^ i becomes diagonal. The 
approximate expressions for 


My and are given in the Appendix IA.2I We only note 


that My 


> My > My 

A2 A3 


with My 


125 GeV, M^^ ~ 50 — 80 GeV and M^^ 


7 keV 


.1 


3 Relic Density Calculation of Two Component Nonther- 
mal Dark Matter 

In the present scenario both the heavier and lighter dark matter components would be produced 
nonthermally in the early Universe. We assume that, when the temperature of the Universe 
was above that of electroweak symmetry breaking, there was no source of production of the DM 
particles and After the breaking of SU( 2 )lX U(1)y symmetry, the self-annihilation of SM 
particles such as W, Z, Higgs, t-quark and/or decays of the Higgs boson act as primary sources of 
the dark matter particles Xj and Xg • The Feynman diagrams of the above mentioned processes, 
which are relevant for the evolution of the number densities of X 2 and Xg, are shown in Fig. [T] 
We compute the number densities of both the dark matter components X 2 and Xg at the present 
temperature (Tq ~ 10“^^ GeV) by numerically solving two coupled Boltzmann equations which 

®The light mass (A(keV) for the dark matter component Xg arises from an order MeV value of u {cf. Sections 
I113 so long as the parameter Kg lies in the range ^ 2 x 10 ^ to 4 x 10 
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Figure 1: Feynman diagrams for dominant production channels of both the dark matter compo¬ 
nents and 


are given below. 


dz 


2Mpi z \/g^{T) 


^g^{T) 


V ^ /r \(Y — ^ ^ 

1.66M2 g,{T) y Y xj 45 i_gg ^2 


X 


E 




,) iyx'-yy')+i^y 


. x=W,Z,f,H 


\'2 X 2 ^X'g Xg ) j I 


( 13 ) 







dz 


2A-4, ' (y ye,\ ^ 

1.66M2^ ^,(T) Xi^XaXa/l^ X 3 xj ^ 45 i_gg ^2 



( 14 ) 


In Eqs. ffT^ and ffT4]) . (Y^ = -d^) is the comoving nnmber density of (Xg)) 

z and T is the photon temperatnre while s is the entropy density of the Universe 1_|. Here 

the nnmber densities of > X 3 Eire denoted by respectively. Fnrther, Mpi is the Planck 

mass and is given by 


VsJZ) 


g.{T) / l dlng,{T) \ 


(15) 


In the above, gp(T) and gs(T) are the effective degrees of freedom related to the energy density p 
and the entropy density s respectively of the Universe throngh the relations p = 5 fp(T)|^T^, s = 
9 s{T)^T^- Thus Pi, is a function of the stated effective degrees of freedom. The thermal 
averages of decay widths (T) and annihilation cross sections times relative velocities (crv) for 
various processes, that occur in Eqs. [I3] and [TU can be expressed as 


(r 


:) = r> 


Kiiz) 


Xi^X^.XU K 2 {z)' 

{(Z'Vxx^X, X, ) ~ 


J = 2, 3, 




cr. 


xx^X 


xAs- (is , 


j = 2,3, X = W^, Z, f, X,, ■ 


(16) 


(17) 


In Eqs. flT6|) . flT7|) Ki is the modihed Bessel function of order i and s is the Mandelstam 
variable. The decay widths T^^^^ ^ and annihilation cross sections axx^x x, d ~ 2, 3, 


®It is to be noted that in the above two Boltzmann equations (Eqs. (IT^ . (1111) 1 we have neglected a term 
involving (crv^ )Y^ . At an earlier epoch during the initial stage of production of (mainly from the 

decay of ), the number density of was very low and hence this term could be neglected. On the other 
hand, at a later epoch when the Universe cools down to a temperature lower than the mass of , the process 
X 3 X 3 X 2 X 2 "'xill not have any significant contribution to the term (crv^ ^ )Y^ even though the number 

density for y is higher. 
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X = Z, /, Xjj Xj) of fho processes mentioned in the subscripts of T and a are given below: 


Xi^X^.X^. 


9 

X-, XX. 
1 J J 


4M2 


1 - 


Ml ’ 


(18) 




27rs 

6 a 


\ 


S-4M2 

- -{ 9 ^ +T 

s-AMl I ^ 1^1 XX' 


9 

•-'v V V •-'v 




g g is — ) 

X, X, X - X . V y, / 




(19) 


a 


^2 ^2 “*‘^3 ^3 


271S 

2 a 


\ 


S-4M2 

X' 

s - 4M2 


2 

•-'T 


Xo Xo Xo Xc 


+ T 


9 9 

X2 X2 X-^ X^ X3 Xg 


^2 ^2 ^3 ^3 ^2 ^ 2^1 ^ 1^3 ^3 


«<.x,v,(s-XC ) 


(s-Af=_)2+(U_A4_)= 


( 20 ) 


(T. 


ww^x ■ X ■ 
3 3 


9^ 9^ 

•Jyvwx-^ “'Xj Xj Xj 


Is-AMl ( 3 -^ + ^ 


M, 


w 


72713 \ s- 4M^ (s - M2 )2 


( 21 ) 


(T. 


.ZZ-^X . X . 
3 3 


g"^ a^ 

i'zzxj x^. x^. 


- "X (3 - 5| + IBf) 


ISvrs \ s ~ 4M| (s — M2 )2 


( 22 ) 


(J, 


//^X . X . 
3 3 


n, 




IbvTS 


(s - A/J )2 


(23) 


In the above equations gijk and gijki are couplings of the vertices involving the fields i, j, k as 
well as i, j, k, I respectively while Mw, Mz and Mf are the masses of W boson, Z boson and 
fermion / (/ is any SM fermion). All the couplings which are necessary to calculate the decay 
widths and annihilation cross sections are given in the Appendix lA.ll 

Finally, the total relic density (hlTh^) of dark matter in the Universe is given in terms of the 
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normalised Hubble constant h = loo^n^/Mpc/s 

VL^h'^ = . (24) 

Once we obtain the comoving number densities {To), (Tq) of both the dark matter compo¬ 
nents Xj and Xg at the present temperature Tq by numerically solving the two coupled Boltzmann 
equations (Eqs. [l3l[Tl]), the individual relic densities and of each of the compo¬ 

nents can be obtained from mm 

f My \ 

= 2.755 X 10^ f ^ j (To) (^ = 2,3). (25) 

While solving these two coupled equations (Eqs. flT^ . ffT4)) h we have adopted the following 
boundary condition: at the electroweak phase transition temperature, which corresponds to 
z (= ^T^) — 0.83, the number densities of both the dark matter candidates are zero. The plots 
(a—c) in Fig. [2] show the variation of the relic densities of both the dark matter candidates Xj 
and Xg with z (inverse of temperature) for different values of model parameters, namely A 12 , 
A 13 . It is to be noted that these parameters turn out to be of order ~ on account of 

the smallness of the comoving number density required to satisfy the observed dark matter 
relic density in the Universe. All these plots (a—c) of Fig. [2] are drawn for the case with 
= 70 GeV and =7.1 keV. The red solid line in each plot of Fig [2] represents the relic 
density of the heavier dark matter component (x^) while the green and blue solid lines denote 
the relic densities of the lighter DM candidate Xg and the total density of both the dark matter 
components respectively. In plot (a) of Fig. |2]we have chosen the values of A 12 , A 13 in such a way 
that the DM particle Xg becomes the dominant component within the dark sector in terms of 
its contributions towards the total dark matter relic density (Dxh^). However, for the other two 
plots, different sets of chosen values of A 12 and A 13 result in the situations where h? < h? 
(plot b) and ~ ^Xg^^ (plot c). It appears from each plot of Fig. [2] that the relic density 

of each dark matter candidate starts growing from an initial value zero (due to the adopted 
boundary condition discussed earlier), thereafter, as the temperature of the Universe decreases 
{z increases), the relic densities of both the DM components increase since more and more dark 
matter particles are produced by the decay and/or self annihilation of the SM particles. Finally, 
the relic densities of both the dark matter particles saturate to the respective particular values 
at z ~ 10 (corresponding to a temperature T ~ 12 GeV of the Universe) which depend upon 
the values of the parameters A 12 and A 13 . It is also seen from Fig. [2] that in all three cases the 
saturation values of the total relic density of the two dark matter candidates always he within 
the range 0.1172 < < 0.1226, as predicted by the PLANGK experiment [2] at a 68 % G.L. 

In Fig. [3] we show the variation of the relic density of the heavier dark matter component Xg 
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(a) 


(^) ^ ^x. 



(*") ^Xj — ^Xg 

Figure 2: Variation of relic densities of both the dark matter candidates with 2 ;. 

with z for three different values of My.^ = 65, 70, 80 GeV with A 12 = 1.10 x 10“^h 

In the three plots (a—c) of Fig. 01 we show the allowed regions in the parameter space A 12 -A 13 , 
^i 2 "^i 3 and ^ 12-^23 respectively. These plots (a—c of Fig. iD are drawn for 50 GeV < < 

80 GeV, = 7.1 keV, 2 MeV < n < 10 MeV, and lO"® GeV^ < a < 10"^ GeVT The choice 
of the numerical ranges of these parameters will be justified in Sections 01 [ 6 ] and [71 In order to 
avoid the late time decay of the heavier dark matter particle into two lighter DM component 
Xg, we adopt the value of the parameter A 23 < 10“®. The constraint for all the plots of Fig. 01 is 
that the calculations of relic densities with the allowed values of the parameters must satisfy the 
PLANGK result for the total relic density of dark matter in the Universe. It is also to be noted 
that, for plot (a) in Fig. 01 we obtain two allowed regions in the parameter space A 12 -A 13 . The 
nonthermal production of the present dark matter candidates proceeds through two processes. 
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Figure 3: Variation of relic density of the dark matter component with z for three different 
values of = 65, 70, 80 GeV. 


namely the decay of the SM Higgs and pair annihilation of SM particles. But, if the mass of 
the heavier dark matter is higher than M^^/ 2 , then it will not be produced through the decay 
of the SM Higgs boson. In this case the heavier dark matter can only be produced by the pair 
annihilation of SM particles. The lower allowed region in plot (a) between A 12 and A 13 is for the 
case when the mass of the heavier dark matter is less than M ^^/2 while the other region is for 
the case when the nonthermal production of the same is through the pair annihilation of SM 
particles. Needless to mention here that the lighter component can always be produced from 
the decay of the SM Higgs boson. The allowed ranges for the model parameters A 12 , A 13 , 6 * 12 , 


Ai2 

Ai3 

012 

(rad) 

013 

(rad) 

023 

(rad) 

~ (0.18 - 1 . 6 ) 
xlO-^i 

~ (0.56-2.6) 
xlO-9 

~ (0.1 - 2.4) 

X10-23 

~ ( 0.2 - 2 . 0 ) 
xlO -11 

~ (0.5 -8.0) 

xlO-13 


Table 1: Allowed ranges for the model parameters A12, A13, 6*12, 6*13, 6*23. 

^ 13 ) ^ 23 , obtained from the plots a—c of Fig. HI are given in a tabular form (Table [T]). Thus 6 * 13 , 
023 and 012 turn out to be very small angles in the range ~ 10 “^^ rad, 10 “^^ rad and 10 “^^ rad 
respectively, as the off-diagonal elements of the mass squared matrix of Eq. (IT^ are constrained 
to be very small. 


13 

































8.0x10'’^ 

6.0x10'’^ 

S 

a 

4.0x10'’^ 

2 .0x10'’^ 

S-OxlO'^"^ 1.0x10'^^ 1.5x10'^^ 2.0x10'^^ 2.5x10'^^ 

012 (rad) 

(c) 

Figure 4: Allowed ranges of model parameters namely, A 12 , A 13 , 6 * 12 , O 23 and 013. 

4 Domain wall formation from restoration of Z 2 

In this present model, though the residual discrete symmetry is spontaneously broken by 
the VEV u of the scalar held S' 3 , there still remains the possibility that the symmetry could be 
restored at a high temperature. This would result in the formation of domain walls. Since the 
scalar held Si is in thermal equilibrium with the Universe, the interaction between the scalars 
Si and S 3 results in a temperature dependent mass term for the held S 3 . The expression for the 
temperature dependent mass term is given by [ST] . 
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We are considering the epoch of the Universe when T -C Ms^- In this limit, we can write the 
above eqnation (Eq. fl26|) i in the following approximate form 




/Ai3M5,T\ 

V 27r2 ) 



T<^Mh 




^13 y2 
2\/27r2 



MsJT 


3 T 

1 + X- 


8 M; 


Si 


15 T2 
128 



(27) 


From Eq. (jl]) we see that there is a wrong sign mass term —'^u^ of the field S' 3 . Combining this 
bare mass term with the temperatnre dependent mass term given in Eq. fl26ll . we can dehne a 
qnantity which is 

MliT) = (28) 

The discrete symmetry Zg will be restored again if the qnantity Mg^{T) changes its sign and 



^13 


Fignre 5: Variation of discrete symmetry (Z 2 ) restoration temperatnre with the parameter A 13 . 


becomes positive. This will happen only if the qnartic conpling A 13 between the helds Si and 
S 3 (see Eq. (jl])) is positive. For negative valnes of A 13 , which are eqnally allowed, there is 
no possibility that the spontaneonsly broken discrete symmetry Z 2 conld be restored again. 
The temperatnre at which this phase transition occnrs is defined as the symmetry restoration 
temperatnre Tr. In Fig. Hj we show the variation of Tr with the allowed valnes of the conpling 
parameter A 13 . It is seen from Fig. [5]that for the allowed valnes of A 13 , Tr varies from ~ 14 GeV 
to ~ 16.5 GeV. 

We note that onr Tr is qnite close to the discrete symmetry restoration temperatnre ~ 15.2 
GeV obtained in the model of Babn and Mahapatra isa. Therefore, the properties of any domain 
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wall that is formed here will be similar to theirs. In consequence, their discussion carries over 
to our case mutatis mutandis and the energy density of such a domain wall would be too little 
either to overdose the Universe or affect the observed near-isotropy of the CMB so long as the 
VEV u of the scalar held S 3 is bounded from above by ~ 10 MeV. 

5 Direct detection of the dark matter candidates and 

u- 

In the present two component dark matter model both the dark matter particles, namely and 
Xg, scatter off the detector nuclei placed at various underground laboratories. These scattering 
processes occur mainly through the exchange of the SM-like Higgs boson The Feynman 
diagrams for the spin independent scattering of both the dark matter particles with the detector 
nucleon [N) are shown in Fig. [HI The expression for spin independent elastic scattering cross 




Figure 6 : Feynman diagrams for the spin independent elastic scattering between each of the dark 
matter particles and the nucleon {N). 


section between the dark matter component x. (* = 2, 3) and the nucleon N through the exchange 
of Xj is given by 


X N^x^N 
(^SI 


y 

TT \ V J 


\ 



(29) 


where is the reduced mass between the dark matter component x. cind the nucleon N of 
mass . In Fq. fl2^ the coupling term among the helds Xj X^ is represented by the quantity 
S'xjX X while / ~ 0.3 |82] is the usual nucleonic matrix element. The expressions of the coupling 
term S'x^x x for i = 2, 3 are given in Appendix lA.ll 

We have already shown that in the present scenario the dark sector is composed of two 
different scalar helds Xj X 3 whose masses as well as interaction strengths with the SM 
particles are different. Therefore in order to compare the spin independent scattering cross 
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section between each of the dark matter particles and the nucleon, computed using the present 
formalism with the results obtained from the present ongoing direct detection experiments such 
as XENON 100 [TB], LUX [77], one needs to multiply the scattering cross section ag} * by 
the factor ^— HT], where Uy. is the number density of the dark matter candidate 

^X2 +’^X3 ' - ’ 

X ■ ^—^X ■ ^ 

X. ■ The rescaling of agi " is due to the fact that the exclusion plots obtained from various 
dark matter direct detection experiments are computed by assuming that the Universe contains 
only one type of dark matter particle, which is not a valid assumption for the present scenario. 
The spin independent scattering cross sections of both the dark matter particles Xj) X 3 given 
in Table [2l It is seen from Table [2] that the couplings for the vertices X^ X^ X^ X^ 


Dark matter 

Mass 




candidate (yj 

(GeV) 

(GeV) 

(pb) 


T 2 

50-80 

(0.44 - 3.9) X 10-9 

(0.02 - 4.42) X 10-27 

10-9 - 10 -® 

X 3 

7.1 X 10-® 

(1.4-6.4) X 10-^ 

(0.14-3.03) X 10-^9 

~ 1 


Table 2: Spin independent scattering cross sections of the dark matter particles and 

are extremely small (due to the nonthermal origin of x^ cind x^)- Hence the spin independent 
scattering cross sections of both the dark matter particles and lie well below the present 
limit which is ~ 10 “® pb to 10 “^° pb for a dark matter particle having mass in the range 
10 GeV- 100 GeV [77]. 

6 7 -ray excess at 1-3 GeV energies from Galactic Centre. 

Our endeavour in this section is to explain, within the framework of our proposed two component 
dark matter model, the 7 -ray excess from the Galactic Gentre observed at an energy range 1 
to 3 GeV by the Fermi-Large Area Telescope. One of our dark matter components, namely y^ 
having a mass in the range ~ 50—80 GeV, decays predominantly into the bb hnal state since 
it has nonzero mixing with the SM-like Higgs boson y . The y bb coupling (g , / = b) can 
be read off from the expression given in Appendix lA.ll It is noteworthy that the small strength 
of the y^y^y^ coupling makes the pair annihilation of two y^’s into bb (through an s channel 
y^ exchange) a negligible competitor. We believe that ours is the hrst model explaining the 
Fermi-LAT 7 -ray excess purely from the decay y^ —)■ bb. The b-quarks resulting from the decay 
of y^ hadronise to produce 7 -rays, the spectrum of which should explain the 1-3 GeV excess 
as seen by Fermi-LAT. The differential 7 -ray flux due to the decay of the component y^ in the 
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( 30 ) 


channel Xj bb in onr two component DM model is given by [83] 


dVLdE 


re PQ 

Att Mv 


J r' 


dN'^ 


dE ■ 


dN° 


In Eq. flHnjl . is tPe energy spectrnm of the photons prodnced with energy E from the 

" ,dN^ 


hadronisation of b qnarkqj. We have nsed the nnmerical valnes of the photon spectrnm (-^) for 
different photon energies given in Ref. [83]. The valne of the dark matter density at the solar 
location, namely Pq, is taken to be 0.3 GeV/cm^ while r© ~ 8.5 Kpc is the distance between the 
GC and the solar location. We have averaged the J factor over the opening solid angles: 

4 


J = 


AQ 



db dl cosb J{1, b) 


with 


J{l,b) = / - 


ds /p(r) 

l.o.s V Po 


and 


= A dl db cos b , 


(31) 


(32) 


(33) 


r = (rQ + — 2 TqS cos6cos/) 


1/2 


(34) 


In Eqs. flTT]) . fl55]) . flMjl . I and b denote the galactic longitnde and latitnde respectively and s is 
the line of sight distance. While compnting the valnes of J, we have performed the integral over 
a region which is situated within an angular distance of 5”|I21 around the GG. The integral over 
s in Eq. fl32p is along the line of sight (l.o.s). In the expression for differential y-ray flux (Eq. 
(jSni)) the quantity T'^ ^^6 is the product of the decay width of the channel —)■ bb and the 
contribution of the component y to the total dark matter relic density (Gt/i^): 


T' 


X,^bb ; 


where 


^Xn ~ 






(35) 


(36) 


is the fractional relic density for the component y . The use of the modihed decay width T' .uu 

2 _ ^2 

(Eq. [35|), instead of the actual decay width of the channel —)■ bb (T^^ .bb) in Eq. (EOj), 


^These originate from the decay of the dark matter component y through the process y —> bb. 
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is required here as we are dealing with a two component dark matter model and the fractional 
amount of the relevant component should be considered for the computation. Needless to mention 
here that, if the entire dark sector consists of only one type of particle (say Xj); then = 1 , 
and consequently V .uc and F .uc are identical. We calculate the decay width of the DM 

A2 _ A2 

component for the bb final state (see Fig. [7] for the Feynman diagram of this decay process) 
and the expression for this decay process is given below: 


X2 


. 



Figure 7: Feynman diagram for the decay channel —)■ bb. 



^•bb — 


n G 

-^^=^(sin 612 cos 023 + cos 0 i 2 sin 023 sin , 

4v27r ^ 


(37) 


where 



In the above, Gp = with v as dehned in Eq. ([T]), is the Fermi constant and is the mass 
of the b quark. 

For the computation of the 7 -ray flux using Eq. flHOD . the astrophysical input is the variation 
of dark matter density (as a function of the radial distance r) in the neighbourhood of the GC. 
This functional relation between p{r) and r is known as the halo profile of DM. In this work, the 
computation of the gamma flux is done considering the NFW profile ca. The general expression 
of the NFW prohle is given by. 


Pnfw 



where 7 is a parameter (index) of order one. In the present calculation 7 
the halo prohle reduces to the form 


Pnfw 





(38) 


1 is adopted. Now 


(39) 
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where = 20 Kpc. The normalisation constant ps is obtained by demanding that the dark 
matter density at the solar location (r = r©) is 0.3 GeV/cm^. 




Figure 8: Left Panel: y-ray flux for three different values of = 50, 60, 70 GeV. Right Panel: 
y-ray flux obtained from the decay of a 70 GeV dark matter particle (y^) for three different 
values of modihed decay width of the channel ^2 bb. 

The gamma-ray flux is calculated following the above discussion and the results are shown 
in Fig. [Hi The Fermi-LAT observational results (black vertical lines) are shown in Fig. jS] for 
comparison. In the left panel of Fig. [HI we plot the y flux from the Galactic Gentre region for 
three values of namely = 50, 60, 70 GeV. From the left panel of Fig. [Hi it is evident that 
the Fermi-LAT results are best described for the choice = 70 GeV. In the right panel of Fig. 
[Hi we adopt the value = 70 GeV and compare our calculated results (for the y flux) with 
the Femi-LAT observational data with three different values of modified decay widths F'^ = 

3.90x10“^^ s“^ (blue dashed line), 4.15x10“^^ s“^ (red solid line) and 4.30x10“^^ s“^ (green 
dashed line) respectively for the process ^2 bb. The right panel of Fig. [HI shows that the 
decay width of F'^ = 4.15 x 10“^^ s“^ best represents the observational results from Fermi- 

LAT. Here we note that obtaining F' .uc in the right ballpark of 3.90x10“^^ s“^ to 4.30x10“^^ 

A2 

s“b oiie requires the soft breaking parameter a to be in the range 10“® GeV^ a ^ 10“^ GeV^. 

We had mentioned earlier that we have taken the canonical NFW profile where the index y 
equals 1 in Eq. fl39p . But in Ref. [T2|, where a model independent analysis of Fermi-LAT data 
with dark matter annihilation has been made to explain the gamma excess in the energy region 
1-3 GeV, a prohle steeper than the canonical NFW prohle with y = 1.26 was taken. The index 
y ~ 1.2 was also required for explaining, within the framework of dark matter annihilation, the 
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Figure 9: 7 -ray flux obtained from the decay of a 70 GeV dark matter particle (y^) for a steeper 
NFW profile with 7 = 1.26. Three lines represent the 7 -ray flux for three different values of 
modified decay width of the channel —)■ bb. 

gamma signal from the Fermi bubble region in Ref. HD]. We therefore also compute the gamma 
flux from the Galactic Gentre in the present framework using a steeper NFW profile with index 
7 = 1.26. The results are shown in Fig. [9] where we plot the calculated gamma ray flux (taking 
an NFW type profile where 7 = 1.26) with = 70 GeV and three values of the modified decay 
width F' .7. In this case the Fermi-LAT data are best represented when F' ,7 = 2.54 x 10“^^ 


7 3.55 keV X-ray line 

We mentioned earlier that our present two component dark matter model contains a scalar 
particle y^, with a mass of order keV, which possesses a tiny mixing with the SM-like Higgs boson 
y^. The two photon decay mode of the dark matter candidate y^ produces monoenergetic photons 
(X-rays) which have been detected by the X-ray telescopes of the XMM-Newton observatory. It 
has been reported in Refs. P in IBB] that, in order to produce the observed X-ray flux from the 
decay of a dark matter particle, the corresponding decay width for the channel DM —)■ 77 must 
be in the range ~ 2.5x10“^® s“^ to 2.5 xl0“^® s“^. In the present scenario, the above constraint 
should be applied on the modified decay width instead of the actual decay width (F;;^^^.^.^) for the 
channel y^ 77 , since we are working in a framework where the entire dark sector is composed 
of two different scalar fields y and y . The modified decay width (F' ,_) of y for the channel 
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Xg —^ 77 is defined as, 


r' 

^ X,^77 



X3^'77 


(40) 


with 



(41) 


as the fractional contribution of the DM component Xg to the total relic density The 

Feynman diagrams for the decay channel Xg 77 ai'e shown in Fig. [101 The expression for the 



Figure 10: One loop Feynman diagrams for the decay channel Xg 77 - 


decay width for the channel x ^ 77 , which takes place at one loop, is given by 


^3 

Gf 


X,^77 


128x/2 


TT'" 


(sin 6^12 sin6'23 - cos 012 cos 023 sin 0 i 3 )^|Fp, 


where 


(42) 


F = F„(T„) + Y,n,,Q}F,(T,) (43) 

/ 

with 

4M^ _ 4M| 

^3 ^3 

Fh/(v) = 2 + 3v+3v(2-v)g(v), 

^/(^/) = -2^/[l + (l-^/)0(^/)]> 

g(r) = arcsin^[r“^^^] . (44) 


Here ttem — 157 is the hne structure constant while Qf and ricf are the electric charge and the 
colour charge of the fermion (/) § involved in the fermionic loop of Fig. [TUI 

®/ is any SM fermion 
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Figure 11: Left panel: Variation of the modified decay width of the decay channel 

X 3 —>■ 77 with A 13 . Right panel: Allowed region in the n —A 13 plane which produces the observed 
X-ray flux as well as satishes the PLANCK limit on the total relic density. 


In the left panel of Fig. [TT]we show the allowed range 1.2 x 10 “® ^ A 13 ^ 2.6 x 10 “® of the 
parameter A 13 for which the modihed decay width (F'^^^.^.^) for the decay channel —)■ 77 lies 
in the range 2.5x10“^® < < 2.5 xl0“^® which is necessary to produce the observed 

X-ray flux from the extragalactic sources such as Perseus, Andromeda etc. The right panel of 
Fig. [m shows the allowed region in the u — A 13 plane. From this plot (right panel of Fig. [TT |1 
one notices that, in order to produce the observed X-ray flux, the VEV u of the scalar held S '3 
must be > 2.4 MeV. The upper bound {u < 10 MeV) on the allowed values of u comes from the 
domain wall constraint (see Section 3.4 of Ref. ISH) which arises from the spontaneous breaking 
of the discrete symmetry Z 2 . Needless to mention here that, for all the points in both panels of 
Fig. El the PLANCK limit on total relic density of the dark matter candidates is always 
satished. 


8 Conclusion 

In this paper we have presented a two-component model of nonthermal dark matter by postulat¬ 
ing two additional SM-singlet scalar helds S 2 and S '3 which interact among themselves and with 
the SM Higgs doublet H. These interaction terms have been restricted by assuming a suitable 
Z 2 X Z 2 discrete symmetry which is softly and explicitly broken to a residual Z 2 symmetry. The 
latter gets spontaneously broken when S '3 develops a VEV of order MeV. Our physical scalar 
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spectrum comprises three particles: (1) the SM-like Higgs boson with ~ 125 GeV, (2) a 

moderately heavy scalar dark matter particle with 50 GeV < < 80 GeV and (3) a very 

light scalar DM particle with ~ 7 keV. We have discussed the issue of domain wall for¬ 
mation which arises in our model due to spontaneous breaking of the discrete symmetry Z 2 and 
have explained why it is unimportant for us. We have computed the 7 -ray flux from the decay 
channel y^ —)■ bb of the heavier dark matter candidate y^ while the X-ray line is generated from 
the decay mode of the lighter DM candidate y^ into a pair of two keV energy photons. These two 
decay channels (y^ —)■ bb, y^ —)■ 77 ) exist due to the fact that both the DM candidates in the 
present scenario possess tiny amounts of mixing with the SM-like Higgs boson y^. We have found 
that the 7 -ray flux originating from the decay of a 70 GeV dark matter particle (y^) into a bb 
hnal state at the Galactic Gentre with a modihed decay width T'^ = 4.15 x 10“^^ s“^ hts the 

Fermi-LAT data well. Finally, we have shown that the modihed decay width of the lighter dark 
matter component y^ for the channel y^ —)■ 77 lies in the appropriate range of 2.5 x 10“^® s“^ 
to 2.5 X 10“^® s“^ as long as the VEV u is bounded from below by 2.4 MeV. 
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A Appendix 

A.l Couplings of the physical scalars y , and y^ 

The couplings of the physical scalars > T3 ~ among themselves and with other SM particles 

are — given below in the limit when all three mixing angles ^12, ^13, ^23 are extremely small. 

— ~^23(3^126 'i3^i ~ 36 'i26'i3Ai2 — 20i26'i3Ai3 + 46'i26'i3A23)- — , 

^XjX^XgXg ~ ~^23(~3012^13 /xi + 30i2^13Ai2 + 20i2013Ai3 — 4012013A23)- — , 

^'xjXjXgXg — “-^23 — 6'23(36*i26*i3fi;2 — 36*i26'i3Ai2 + 36*126*13 Ai3 — 66'i26'i3A23) , 

9x,x,x, - - ^13mAi3 , 

^XjXjXj ~ ^12'i'Ai2 — 6'23(2 mA 23 — 6*13'l'Ai2) , 

^'xgXgXg ” “ ^3 T ^ 13A13 “ 6^23(“26^ 1 26 ^ 13MA13 T 6^ 12UA13 T 46*i26^ 13UA23) , 

5x^X2 Xj ~ “^23 (66*126*13^/1;^ — 46*126*13 uAi 2 — 2012mAi 3 — 26*i26*i3uAi3 + 46*12 mA 23) 

— uAi2 — 26 *i 3 mA 23 

^'x^XgXg — “^23 ( —60126 *i3U/xj + 46*i26*i3uAi2 + 26 *i2mAi3 + 26*i26*i3uAi3 — 46*i2uA23) 

— 'i^Ai3 + 20 i3uAi3 — 39 i 3UK,^ 

9x^x^x^ — ~^23 (3 m/x 3 + 20 i 3 uAi 2 — 36*i3uAi3 — 4MA23) ~ 26 *i 26 *i 3 uAi 3 + 6 *i 2 uAi 3 

+ 46*i26*i3mA23 , 

Qx^x^x^ - “^23 {QOi^UK^ - 2v\i2 - 46*i3mAi 3 + 2uAi3 - 46 *i3mA23) - 66 *i 26 *i 3 U/x^ 

+ 40 i 2^13 uAi 2 + 20 i 2 'uAi 3 + 20i2^13uAi3 — 40i2'uA23 , 


9wWx-^ 

rv-* 

2M^ 

V ’ 

9wwx^ 

rv-* 

-^(012 + ^13^23 
V 

9wwx^ 

rv-* 

(n a a 

1(713 (712(723 

V 

9zzx^ 

rv-* 

M^w 

V 

SzZXg 


(012 + 013 023) 
V 

9zzx^ 

rv-* 

-^(013 - 012023) 

V 

9ffx^ 

9ffx2 


Mf 


V 

Mf 

-(012 + 013023) , 

V 

^//Xg 


-^(013 — 012 023) • 

V 



A.2 Masses of physical scalars y , and y^ 

The masses of the physical scalar and Xg in the limit of extremely small 6 ^ 12 , 6 ^ 13 , 623 are 

given as fellows: 

= (2k^v‘^ + AuvOisXis + 2 q ; 0 i 20 i 3 ) ^ + £^(^ 12 ) + ^^(^ 23 ) + ^(^ 13 ) ) 

= {Pj + Ai2n^ + 2A23'U^ + 2023 (o^ ~ P2^126'l3 + 2n^Kg0i2013 ~ n^Ai2012013 
— 2nnAi30i2 — 2n^A230i20i3) } ^ + ^(^ 12 ) + ^(^ 23 ) + ^(^ 13 ) 5 

M^3 = 12^3^^ — 2q;0i 2013 — 4Mn0i3Ai3 — 2023 (ct ~ P2^126 'i 3 + 2n^Kj012013 

—n^Ai2012013 ~ 2MnAi30i2 — 2 m^A23012013) } ^ + ^(^ 12 ) + ^(^ 23 ) + ^(^ 13 ) • 
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